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ABSTRACT 



We have been monitoring the supernova remnant (SNR) 1987A with Chandra 
observations since 1999. Here we report on the latest change in the soft X-ray 
light curve of SNR 1987A. For the last ~1.5 yr (since day ~8000), the soft X- 
ray flux has significantly flattened, staying (within uncertainties) at fx ~ 5.7 x 
10~ 12 erg cm~ 2 s _1 (corresponding to Lx ~ 3.6 x 10 36 erg s _1 ) in the 0.5-2 keV 
band. This remarkable change in the recent soft X-ray light curve suggests that 
the forward shock is now interacting with a decreasing density structure, after 
interacting with an increasing density gradient over ~10 yr prior to day ~8000. 
Possibilities may include the case that the shock is now propagating beyond a 
density peak of the inner ring. We briefly discuss some possible implications on 
the nature of the progenitor and the future prospects of our Chandra monitoring 
observations. 

Subject headings: ISM: supernova remnants — supernovae: individual: SN 1987A - 
X-rays: individual: SN 1987A 
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1. INTRODUCTION 

The known distance (d ~ 50 kpc to the Large Magellanic Cloud [LMC]) and 
the identification of the progenitor (Sk -69°202, a B3 I star) make SN 1987A a 
uniquely useful object for a detailed study of core-collapse supernova and the earliest 
stages of evolution of a supernova remnant (SNR). For the last ~10 yr, SNR 1987A 
has been dominated by emission from the blast wave interacting with the dense 
circumstellar medium (CSM), dubbed "the inner ring", which was produced by 
interaction between the slow equatorial wind from the progenitor's red supergiant (RSG) 
phase and the fast spherical wind from the subsequent blue supergiant (BSG) phase 



(IBurrows et al. 1995 



Lundqvist fc Fransson 1991 



Luo fc McCrav 19911 ). We have been 



monitoring SNR 1987A with Chandra since 1999. Thanks to the high angular resolution, 
our Chandra monitoring observations have allowed us to reveal unprecedented details of the 
dynamical and thermal evolution of the hot gas in SNR 1987A (e.g., Burrows et al. 2000; 
Michael et al. 2002; Park et al. 2002,2004,2005,2006; Racusin et al. 2009). SNR 1987A has 
also been observed by XMM-Newton (e.g., Haberl et al. 2006; Sturm et al. 2010). 

The evolution of SNR 1987A has been highlighted by two characteristic events: (1) 
the emergence of optically-bright hot spots since 1997 (day ~ 3700) and (2) the up-turn in 
the soft X-ray light curve that was coincident with a deceleration in the radial expansion 
rate in 2004 (day ~ 6000). The former indicated the start of the interaction between 
the blast w ave and the dense p rotrusions that had grown inward along the surface of the 



inner ring (IMichael et al. 20001 ) . The latt er was interprete d as the event of the blast wave 



reaching the main body of the inner ring (IPark et al. 20051 ). The continuous propagation of 



the blast wave into and/or beyond the inner ring would provide an excellent observational 
opportunity to study chemical and density structures of the progenitor's stellar winds and 
thus to reveal its detailed late-stage evolution. 
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Here we report the latest soft X-ray light curve from our Chandra monitoring 
observations. It is evident that the soft X-ray flux from SNR 1987A has been nearly 
constant over the last ~1.5 yr. This recent X-ray light curve is in contrast to the rapid 
brightening for the prior ~10 yr. This remarkable change in the light curve indicates that 
SNR 1987A may now be entering yet another evolutionary phase. In this Letter, we report 
this new behavior in the soft X-ray light curve and briefly discuss possible implications. 
This work is based on ACIS and HETGS zeroth-order observations that we have performed 
as of 2010 September. The detailed analysis of the first-order HETGS spectrum obtained 
by our monitoring observations will be presented elsewhere. 



2. OBSERVATIONS & DATA REDUCTION 

As of 2010 September, we have performed a total of twenty-two monitoring observations 
of SNR 1987A with the Advanced CCD Imaging Spectrometer (ACIS) and the High Energy 
Transmission Gratings Spectrometer (HETGS) on board Chandra. To mitigate increasing 
photon pileup effects in the ACIS data due to the continuous brightening of SNR 1987A, we 
have switched our monitoring instrument configuration from the ACIS to the HETGS since 

2008 July (day ~ 7800). The soft X-ray flux increase has been ~10— 15% per ~6-month, 
and the cross-calibration of the observed flux between these configurations is critical to 
make a smooth, reliable transition. Thus, we have performed contemporaneous observations 
using both the ACIS and HETGS (separated by <^ one week) for two epochs (2008 July and 

2009 January). During the same time period, the correction tool for the charge transfer 
inefficiency (CTI) in the ACIS data has become available as part of the standard Chandra 
data reduction software CIAO. Thus, we adopted this standard CTI-correction instead 
of our previous method developed by Townsley et al. (2002). We carefully investigated 
the effect caused by these new instrument setup and CTI-correction method using the 
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zeroth-order HETGS and bare ACIS data. The net effect is small (~3% discrepancy in 
the observed 0.5-2 keV band flux, primarily caused by the new CTI-correction). This 
correction factor is accounted for in our updated X-ray light curves. We re-generated 
auxiliary response files for all twenty-two epochs using CALDB version 4.3.0. Other than 
these changes, we used the same data reduction process as that used in our previous work. 

In the spectral analysis to meas ure the X-ray fluxes, w e used the non-equilibrium 
ionization (NEI) plane shock model (IBorkowski et al. 20011 ). We implemented the NEI 



version 2 in XSPEC with an augmented plasma model using the updated atomic data 
that account for the inner-shell lineal . To be consistent with our previous works, we 
used a two-component NEI shock model to fit the observed X-ray spectra. We varied 
the electron temperature, ionization timescale, and normalization for each of the soft 
and hard component in the model fits. We fixed elemental abundances for N, O, Ne, 
Mg, Si, S, and Fe at the best-fit values estimated by our deep Chandra grating data 



f lZhekov et al. 2006 



Zhekov et al. 20091 ) . We fixed the abso rbing column at Nn = 2.35 



x 10 21 cm -2 based on our previous work (IPark et al. 20061 ). Our Chandra monitoring 
observations are summarized in Table [TJ The observed X-ray fluxes for all Chandra 
monitoring observations are presented in Table El We e stimated the ACIS pileup effect 



based on the event grade distribution ( IPark et al. 20071 ). In Table [2J the observed 0.5-2 



keV band fluxes have been corrected for the pileup using the correction factor shown in 
Table HJ 



x The updated model has been provided by K. Borkowski: cf. the discussion by Badenes 
et al. 2006. 
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3. Light Curve and Expansion Rate 

The latest X-ray light curves of SNR 1987A are shown in Figure [Ik. As of 2010 
September, the observed 0.5-2 keV band X-ray flux is fx ~ 57.9 x 10 -13 erg cm -2 s _1 (Lx 
~ 3.63 x 10 36 erg s _1 at d = 50 kpc, after correcting for N-r = 2.35 x 10 21 cm" 2 ), which is 
~40 times brighter than it was in 1999 and nearly three orders of magnitude brighter than 
in 1992. While the X-ray flux has been significantly increasing over the last ~20 yr, the 
latest data show that the X-ray flux increase rate has been significantly reduced. This is 
clearly seen in the soft band (0.5-2 keV) light curve (Figure [Tk.). For the last ~1.5 yr, the 
0.5-2 keV band X-ray flux has been nearly constant at fx ~ 57 x 10~ 13 erg cm -2 s -1 (Lx 
~ 3.6 x 10 36 erg s _1 ). This change is unclear in the 3-10 keV band, probably because of 
a lower flux increase rate and/or a physically different origin than the soft band emission. 
Upcoming monitoring observations are needed to reveal the nature of the hard band light 
curve. Nonetheless, the flattening of the soft X-ray light curve is compelling, and we discuss 
only the soft band light curve in this work. 

Figure [lb shows the updated radial expansion rate. The radius of individual images (in 
the 0.3-8 keV band) and the expansion rate are estimated by the same methods described 
in Racusin et al. (2009). We tested any systematic discrepancy in the radius measurements 
between the ACIS and HETGS zeroth-order images using our cross-calibration data taken 
in 2008 July and 2009 January. We have also compared our radius estimates between 
the previous and current CTI corrections. These calibration results showed negligible 
discrepancies in the radius estimates (<^1%). Thus, for the data taken since 2008 July, 
we present the radial expansion rate based on the HETGS zeroth-order images with the 
standard CTI-correction method. The radial expansion rate has been consistent at v ~ 
1600 km s~ x for the last ~6 yr. 
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Discussion 



Ever since the blast wave started interacting with dense protrusions of the inner ring 
at day ~3700, the soft X-ray light curve showed a rapid flux increase (e.g., Park et al. 
2004). During days ^6000-65 00, the soft X-ray light curve turned up to show a steeper 



flux increase rate than before (IPark et al. 20051 ). This phase was interpreted such that the 



evolution of SIN 



R 1987A entered a stage in which the shock is reaching the main body of 



the inner ring (IPark et al. 20051 ). Such an interpretation was supported by th e significant 
deceleration in the radial expansion rate at day ~6000 (IRacusin et al. 20091 ). Other 
observational supports include d the substantial decrease in the line broadening measurec 



by deep Chandra grating data (IZhekov et al. 2005 



Zhekov et al. 2009 



Dewev et al. 2008h . 



the prevalence of optical spots all aro und the inner ring (IMcCray 20051 ). and the significant 
IR flux increase (IBouchet et al. 20051 ) at the same epoch. These evolutions of the X-ray 
light curve, radial expansion rate, and individual line widths suggested that X-ray emission 
was dominated by the blast wave interacting with the low-density HII region before the 
inner ring until day ~4000, and that since then soft X-rays primarily originate from 
the blast wave interacting with the dense inner ring. Recently, line fluxes from some 
i ndividual optical spots s howed a hint for flattening of the optical light curve at day ~7000 



(jGroningsson et al. 20081 ). 



Our Chandra monitoring observations show that the flux increase rate appears to be 
lower at days ~6500-8000 than at days ~6000-6500 (i.e., fx ~ t 6 as opposed to ~ t 10 for 
days ~6000-6500). Subsequently, the soft X-ray flux has been nearly constant for ~1.5 
yr (Figure [IK). The latest flattening of the X-ray light curve indicates that SNR 1987A 
is now entering another phase. Assuming that the observed X-ray flux is dominated by 
the interacting density, the X-ray flux can be considered to be proportional to the volume 
emission measure EM ~ n\ V, where n e is the electron density and V is the emission 
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volume (Although there was a deceleration in the expansion rate at day ~6000, the shock 
velocity has been constant sin ce then. Also, the tem perature dependence of the X-ray 



flux is small: i.e., /x oc T 0,6 (IMcKee & Cowie 19771 )). The exact geometry of the X-ray 



emitting volume in SNR 1987A is uncertain. Previous works suggested that the bulk of 
X-ray emission originates from a disk-like shell between the forward and reverse shocks, 
which encompasses the inner ring. An upper limit for such an emitting volume may be 
approximated by a spherical shell with an expanding thickness, thus V oc r 3 where r is the 
SNR radius. A lower limit of the X-ray emitting volume could be a circular disk with a 
small thickness, thus V oar. These two cases can serve as conservative limits that are useful 
to approximate the temporal evolution of the X-ray emitting volume. The reality is likely 
in between these limits. For V oc r s where s = 1-3, the X-ray flux would be fx oc EM oc 



nl r s 



. Assuming that n e oc r a , we have fx oc r 2a+s . If we fit the light curve as a power law 



(PL), fx oc t@, then since the expansion is linear (Figure [Tb), fx oc r 13 , and therefore (3 = 
2a + s. This simple relation suggests that, for (3 ~ 6 during day ~ 6500-8000, the blast 
wave was still interacting with a radially-increasing density (n e oc r 1,5 for s = 3, and oc r 2,5 
for s = 1), but at a lower rate than that (n e oc r 3 ' 5 for s = 3, and oc r 4 ' 5 for s = 1) during 
day ~ 6000-6500. Since day ~8000 the flux has been constant, which corresponds to the 
shock interacting with a decreasing radial density profile (n e oc r~ for s — 3, and oc r~ ' 5 
for s — 1). We note that our PL fit analysis is likely an oversimplification. We intend to 
provide insights on physical implications imposed by the recent flattening of the light curve. 
Quantitative interpretations inferred from these simple PL fits before day ~8000 are far 
from conclusive. 

These simple estimates of the density structure into which the blast wave is now 
propagating suggest that a significant fraction of X-ray emission is now coming from the 
shock interacting with a CSM having lower density than the densest region along the inner 
ring. Possible origins for such a change in the ambient density may include an intriguing 
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scenario that the blast wave might have reached a density peak of the inner ring at day 
~8000, and that it is now propagating into a decreasing radial density beyond the inner 
ring. A likely scenario for the progenitor's evolution before the explosion of SN 1987A is a 
BSG phase that had gone through an earlier RSG stage. The standard stellar evolutionary 
models predict an RSG wind with a constant velocity and a radially-decreasing density 



profile (to ~ r 2 ). Our estimates of the latest radial density profile (to ~ r 15 - r~ 



■0.5\ 



suggest that the blast wave of SNR 1987A may be entering the RSG wind beyond the 
interacting layer between the BSG and RSG winds. If this is true, the soft X-ray light curve 
would eventually decrease at a rate of fx oc t^ 1 - t~ 3 (depending on s) as the shock would 
be interacting with the RSG wind (to ~ r~ 2 ). If the overall RSG wind is spherical beyond 
the inner ring, f x oc t^ 1 is expected, which is generally consistent with standard SN models 



flChevalier 19821 ) . 



We note that Borkowski et al. (1997) predicted the soft X-ray light curve of SNR 
1987A based on 2-dimensional hydrodynamic simulations. In their models, the soft X-rays 
are dominated by emission from dense gas that has been heated by shocks transmitted into 
the CSM ring as it is overtaken by the blast wave. The resulting soft X-ray light curves are 
sensitive to the assumed gas density, too, in the inner ring. The X-ray flux in their model 
for no = 32000 amu cm -3 reaches a peak value of fx = 6 x 10~ 12 erg cm" 2 s -1 at t — 17 
yr, which roughly agrees with the observed results shown in Figure [T^, (fx — 5.7 x 10~ 12 
erg cm -2 s" 1 at t = 23 yr). However, the actual hydrodynamics of the interaction are more 
complex than such simple simulations. While Borkowski et al. (1997) assumed that the 
CSM ring was a circular torus of constant density, the inner ring is not a circular torus, 
having many inward protrusions manifested by the optical spots (e.g. , Sugerman et al. 



2002). The density in the inner ring spans a range of a factor of ~30 (iMattila et al. 20101 ). 
Nonetheless, these results suggest that more realistic 3-dimensional simulations for the 
interaction of the blast wave and the inner ring would help to reveal the detailed nature of 



- 10 - 



the X-ray light curve. 

Other models for SN 1987A include cases that th e progenitor might have gone through 



a Luminous Blue Variable (LBV) phase (ISmith 20071 ). which most likely involved episodic 



eruptions of strong asymmetric stellar winds. In this scenario, the inner ring might have 
been the relic dense CSM produced by the latest LBV eruption before the SN. The radial 
density profile beyond the inner ring then might have a periodic "density wave". Such a 
density structure may result in a periodic variation in the soft X-ray light curve in the 
future. 

We also consider the latest soft X-ray light cu rve within the frame work of our reflected 



shock structure (RSS) model (jZhekov et al. 20101 ) . The current radial expansion rate 



appears to be well-fitted by the RSS model, as the expansion rate stays the same until day 
~8600. The RSS model fit shows that the contribution in the soft X-ray emission from 
the transmitted shock into the densest region of the inner ring peaks at day ~6700-7800, 
and then it significantly decreases after day ~8000 (Figure [lb). These results generally 
support that soft X-ray emission originating from a low density region has recently become 
significant. 

Although the scatter in the radius measurements appears to become large for the last 
~2 yr, we do not see clear evidence for an acceleration in the expansion as of 2010 September 
(Figure [lb). This is probably because the shock is in transition between the inner ring and 
the RSG wind, and our data may not yet be sensitive to an accelerated expansion rate. 
Or, the shock may still be interacting with the inner ring (see above discussion). Future 
Chandra observations are essential to monitor any change in the expansion rate as well as 
in X-ray light curves. We note that the X-ray radius of SNR 1987A (~0"78 at day ~8600) 
is slightly smaller than that of the optical inner ring (~0'.'83, Jakobsen et al. [1991]). Our 
radius estimates are based on a simple image model (assuming 4-lobes and an underlying 
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torus) fit to deconvo lved images to measure the radius of the peak intensity of the SNR 



(IRacusin et al. 20091). Other methods have bee n used to meas ure the size of the SNR in 



optical (IJakobsen et al. 19911 ) and radio bands ( INg et al. 20091 ). Thus, it is not surprising 
to see systematic offsets among the radius measurements in different wavelengths. Ng et 
al. (2009) confirmed that there is a systematic discrepancy (by ~20 %) between the radio 
and X-ray measurements of the radius of SNR 1987A because of the different methods used 
for X-ray and radio data. Therefore, direct comparisons of the SNR's angular size between 
different wavelength bands are difficult. 

Our continuing Chandra monitoring observations of the X-ray flux and the expansion 
rate are essential to reveal the true origin for the dynamically changing soft X-ray light 
curve. The high resolution first-order HETGS data of our monitoring observations will 
be useful to reveal the spectral evolution of SNR 1987A. Particularly, our deep HETGS 
observation of SNR 1987A that is scheduled in 2011 March will be critical to study the 
detailed evolution of individual line broadenings and flux ratios, and thus the SNR's 
evolutionary phase. If the contribution from the reverse shock in the observed X-ray 
emission has substantially increased, we may expect a bleaching-out of optical Ly Q and H a 
lines in near future due to photo-ionization (by X-ray photons) of neutral H just ahead of 
the reverse shock as predicted by Smith et al. (2005). 



The authors thank K. Borkowski for providing the augmented NEI shock model. This 
work was supported in part by SAO under Chandra grants GO9-0082X and GO0-11072X. 
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Table 1. Chandra Observations of SNR 1987A. 
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9447 


1.03 


10853, 10854 














10855 














10222, 10926 


2009 July/ 
September 


8202 


HETG (1.1 s) 


58.2 


8505 


1.03 


12125, 12126 


2010 March 


8429 


HETG (1.0 s) 


63.9 


9332 


1.03 



11090 
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Table 1 — Continued 



ObsID 


Date 


Age a 


Instrument 13 


Exposure 
(ks) 


Counts 


Correction 


13131, 11091 


2010 September 


8619 


HETG (1.0 s) 


54.4 


8361 


1.03 


a Days since 


the SN explosion. 













b For HETG, the zeroth-order data are used. The frame-time is shown in parentheses. 



c This is the flux correction factor accounting for the CCD pileup and cross-calibration with 
our previous CTI-correction method (by Townsley et al.) that was adopted before 2008 July. 

d For this epoch, simultaneous observations with the ACIS and HETG were performed to 
calibrate fluxes between two detectors. We here present only HETG data. 
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Table 2. The Observed Chandra X-ray flux of SNR 1987A. 



Age a Flux b Flux b 

(days) (0.5-2 keV) (3-10 keV) 



A r*nn 

4609 


1 r i n a 1 

1.52±0.41 





.64±0. 


1 O 

.13 


A T1 1 

4711 


1 n 1 n r 

1.68±0.25 





f7A 1 n 

.79±0. 


on 

30 


5038 


1 n 1 /i 

2.38±0.14 





.84±0. 


no 

.08 


5176 


/"» 1 n a 

2.68±0.43 


1 


no 1 n 

.02±0. 


.27 


5407 


3.57±0.50 


1 


1 r 1 a 

.15±0. 


OO 

.28 


r - r - "i 

5561 


/i on 1 n or 

4.39±0.35 


1 


or in 

,35±0. 


33 


r* 1 

5791 


r - 717 1 n r* 

5.77±0.35 


1 


,66±0. 


.15 


5980 


1 n a 1 

6.88±0.41 


1 


,82±0, 


.16 


6157 


8.70±0.44 


2. 


a in 

,24±0. 


1 

.18 


6359 


1 1 1 1 a n 

12.11±0.61 


2 


r» 1 in 

.21±0. 


.16 


6533 


16.26±0.81 


2 


f~r i n 

,37±0. 


or* 

.26 


1-7-1 

6716 


1 a at i n on 

19.95±0.80 


2 


or i a 

.85±0. 


.26 


6914 


24.23±0.73 


3 


.54±0. 


.18 


7095 


28.72±0.86 


3 


.98±0. 


.24 


7271 


33.76±0.68 


4 


.52±0. 


32 


7446 


37.91±1.14 


4. 


.46±0. 


.27 


7626 


40.89±1.64 


5 


.36±0. 


.59 


7799 


47.93±1.92 


5 


,55±0. 


.44 


7997 


51.10±2.56 


6 


,28±0. 


.69 


8202 


56.07±2.80 


7. 


10±0. 


.92 


8429 


57.25±2.30 


6 


,66±0. 


53 


8619 


57.89±2.90 


7. 


,88±1. 


.02 
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a Days since the SN explosion. 

b Observed X-ray fluxes are in units 
of 10~ 13 erg cm -2 s _1 . Errors are 
with 90% confidence level, estimated 
by "flux error" command in XSPEC. 
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Fig. 1. — (a) X-ray light curves of SNR 1987A from our Chandra monitoring observations. 
The ROSAT data are taken from Hasinger et al. (1996). The shaded area shows the range 
of the predicted light curve assuming that fx oc t _1 -t -3 starting at day ~8600 (see the text), 
(b) Radial expansion measure from the fits to the quadrilateral lobes plus torus model for 
the broken linear fit with an early fast shock velocity (solid line) transitioning to the later 
slower shock velocity (dashed line), (c) The temporal evolution of the fractional flux from the 



